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1. INTRODUCTION 



Decoherence - the phenomenon whereby a quantum system initially in 
a pure state evolves into a improper mixture due to interaction with its 
environment - has been studied extensively, using both toy modelsp| and 
more complex models representing real physical situationsP). Sources of 
decoherence such as scattering 0, |^ and quantum gravity]^, |^ and both 
internal [^, ^ and external [|1^, |rT| environment are taken into consideration. 
0mnes[O] recently presented a general theory of the effect of decoherence 



that includes both the harmonic model of Caldeira-Legett |]TU| and the ex- 



ternal environment considered by Joos and Zeh|jTT[ and is related to the 
quantum state |0 diffusion model. 

Decoherence is especially interesting as it promise to solve the old problem 
of how to derive the classical behaviour of macroscopic bodies from quan- 
tum principles. In fact, in the case of a macroscopic body, decoherence is 
believed to suppress the off-diagonal elements of the spatial reduced matrix 
of the centre of mass, irrespective of initial conditions. This is equivalent 
to establishing a superselection rule for the position of macroscopic bodies 
and to saying that we cannot experience spatial macrosuperpositions. We 
can thus explain why we do not encounter states of this kind when looking 
at everyday objects. To derive such a result it is not necessary to consider 
an external environment. In fact, a system of many particles like a macro- 
scopic body may be considered as consisting formally of a "collective system" 
described by the collective coordinates, and the system described by its mi- 
croscopic coordinates, which can act as an (internal) environment if the two 
formal systems are coupled. This coupling may be either guaranteed by 
some constraints or caused by an external potential. Of course in the case 
of an isolated macroscopic body, the Hamiltonian consist of the sum of two 
separate Hamiltonians relative to the centre of mass and internal variables 
respectively, there is no coupling between them and so decoherence of the 
collective variable is not possible. However, we wish to signal a possible case 
of decoherence, with suppression of correlations in the total momentum, for 
a system of particles even if isolated, provided that we consider relativistic 
corrections to the non-relativistic Hamiltonian. It would pointed out that de- 



coherence of isolated systems has been discussed before in various contexts: 
HaUiwelipi investigated the variables that will generally become effectively 
classical as the local densities (of number of particles, momentum, energy), 
Calzetta and Hufl^ the decoherence of histories of certain correlation func- 



tions in a field theory. We should also mention that the diagonalization in the 
square of momentum basis has been argued for a system of (non-relativistic) 
particles interacting with the linearized gravitational field ||15||. 

It is know that in classical electrodynamic a system of interacting par- 
ticles may be approximately described to terms of second order in 1/c by 
a Hamiltonian depending only on their positions and velocities, since the 
radiation comes in with the terms of order \jc'. The same may be said 



for particles subject to gravitational interaction |T^. In this approximation 
the total momentum of the system is no longer independent, but is coupled 
to the internal degrees of freedom. So, if we adopt the above mentioned 
Hamiltonian (after due symmetrization) as the Hamiltonian operator within 
non-relativistic quantum mechanics, we may expect to have decoherence for 
the collective variable tracing over the microscopic ones. This may occurs 
even if the particles are not interacting, thanks to the relativistic correction 
to their kinetic energy. The coupling we are talking about commute with 
the non-relativistic free Hamiltonian of the centre of mass (cm.). As we 



learned from studies of wide-open quantum systems p!7[] , the diagonalization 
is to be expected in the basis of any dynamical operator commuting with the 
coupling terms and so your candidate to be a preferred basis is that of total 
momentum. The importance of an interaction commuting with the preferred 
basis was already stressed by Zurek [|T^ ]. 

Since the aim of this note is merely to briefly signal such a possibility, we 
will consider for the sake of simplicity a system of n free spinless particles 
of equal mass m in one dimension, a one dimensional ideal gas in fact. We 
stress that we intend to work within the non-relativistic quantum mechanics 
and consider the correction to the kinetic energy and its consequences as a 
relativistic effect. 

The model certainly is not realistic but we intend to offer the result as a 
further case of decoherence with not usual features. Not usual because the 
system is isolated and also because in our case the reduced statistical oper- 
ator diagonalizes on the basis of the total momentum, which thus acquires 
the status of preferred basis. We believe all this is interesting in itself. We 
do like to stress that the calculations will be carried out without recourse to 



any master equation but by following the Schroedinger evolution of a pure 
state. With the kinetic operator used, the Galilean invariance is destroyed, 
so we have to specify the reference frame we are referring to. A Galilean 
transformation corresponds to a translation of momenta. For the approxi- 
mation we are using to make any sense, the reference frames to take into 
consideration are obviously the ones (if they exist) for which the absolute 
value of momentum for each particle is not significantly likely to be of the 
order of or greater than mc. This requisite can be expressed by applying 
the condition {pf) <^ m^c^ for every i. It will certainly be satisfied if the 
more restrictive condition YH=i{Pi) *^ ttl'^c^ also is. To establish a definite 
reference frame, we will adopt the one in which, for a given wave function, 
the last mentioned expression gives the lowest value. It will be the reference 
frame in which the mean value of total momentum is zero. Even selecting 
the optimal reference frame, however, the validity of our approximation is 
not automatically guaranteed and our result will be only significant if the 
condition on the values of (pj) is satisfied. 



2. THE SYSTEM OF FREE PARTICLES IN ONE DIMENSION 
WITH RELATIVISTIC CORRECTION TO THE KINETIC EN- 
ERGY 

The free Hamiltonian which we will consider for quasi-relativistic particles 
is: 

^ 2m ^ 8m3c2 



Using the momentum basis and with ^p{pi,P2, ■ ■ ■ ,Pn) indicating the initial 
state, we have the state at time t as follows: 

2 4 

\x{t)) = e"**/^E> = l(2^~8m^c^)l V{Pl,P2, ■■■,Pn) 

We now change the variables: pi = P/n+rji, where P is the total momentum. 
The rji are not independent, since they must satisfy the relation J27=iVi = 0- 
We will consider as independent variables P together with the first n — 1 rel- 
ative momenta. The last relative momentum is expressed as ?7„ = — J2^=i Vi- 
The Jacobian determinant of the transformation is equal to 1. We will write 
the integration volume element for the new variables as dV = dP dS, where 



dS = drji dri2 ■ ■ ■ drjn-i. Using the new variables and taking into account the 
relation J2^=i ^i = the expression above now becomes: 

it, P^ P* , y^" ,'^i ^i 3P^ 2 ''., .. v^)] ^,- r. 

v(t)) ^ e '' ^2nm Sn-'m'^c^ ' Z-ti = l'- 2m Sm'-^ c'' 4n-^ m'-^ c'' '» 2nm?c^ ' (P ( X^, T/l , . . . , T7„ ) 

where $(P, 771, 772, .. . ,77„) indicates ip{P/n + rii, P/n + 772, . . . , P/n + rjn) Note 
that here and in the various expressions that follow rjn stands for = — J27=i Vi- 
The projection {P'\x(t)) of the state vector at time t on an eigenstate \P') 
of the total momentum is: 

it r P'2 P'* I y^" (-''? ''I F'^ .. n^ -P'. „ n-'^'|] ^/ ,-■/ n 

g fil2nm 8n?m?c^ ^» = 1^2m 8m?c^ 4n^mW^'i 2nm'-^ c'^ * $ ( P , Tl^ , . . . , 77^ ) . 

By tracing (x(t)|P")(P'|x(t)) upon the relative momenta i.e. integrating over 
dS we get the reduced density matrix element on the total momentum basis: 

P ^ ^^ ^ e fi ' 2nm 64m3c2 ' I pi pll 

The factor Jp/p// indicates the integral 

Jp,p. = I <I>(P', 7/1,7/2,..., Vn) $*(P", ^1, ^2, . . . , r/„) e'T:.^^+^vt) dS 

where w = 3{P'^ - P"'^)/{Ahn^rn?c^) and ^ = (P' - P")/{2hnm^c'^). The 
diagonal matrix elements are constant of motion, but when P' 7^ P" the in- 
tegral in the previous expression may cancel out with time, in which case the 
off-diagonal elements are cancelled out, meaning that we have decoherence. 



3. EVALUATION OF THE REDUCED DENSITY MATRIX EL- 
EMENTS 

In order to understand the time behaviour of the matrix elements we are 
interested in is thus necessary to evaluate the integral Ipipn. To take into 
account the constraint rjn = — 127=1 Vi ^^ multiply the integrand by the delta 
function 6{Y17=iVi) using its Fourier representation (l/27r) j e^^i^^'^^'^dk and 
consider rjn too as a independent integration variable. Writing 

^(P', P"; 771, . . . , 77„) = <I>(P', 771, 772, . . . , r/n) ^(P", 771, 772, . . . , r/n) 



and h{rii) = wrjf + zrjf our integral becomes 

Ip,pu = (l/27r) JdkJ e*^''i+**'^(''i)rfr7i / e''''^^+'^^^'^^Ur]2 . . . 
I e*''"+**'^(''")rfr/„ G{P', P"- r]i, r/2, • • • , Vn) 

Let us start with the integration on rjn. We are interested in the behaviour 
of the off diagonal matrix elements at large time, so we shall use asymptotic 
expansions methods to evaluate the result. More precisely, we will look 
at the method of stationary phase [0. According to which, when t is a 
large quantity, the major contribution to the value of the integral arises 
from the immediate vicinity of the end points and from the vicinity of those 
point at which h{r]n) is stationary, i.e. its first derivative h'{rjn) is equal 
to zero. In the first approximation the contribution of stationary points 
is more important. If we suppose that the function representing the intial 
state is square integrable and thus going to zero at infinity, we can disregard 
the end points contribution at all. Since only the neighbourhood matters, 
asymptotic evaluation consists of replacing h{rin) near any stationary point 
r]* as hij]'^) + {l/2)h"{rj'^{rjn — Vn)^- By the same argument, the factor of 
^ith{ri„) jj^ ^YiQ integrand, assuming it is a continuous function, is replaced by 
its value at r^* and the integration is safely extended from — cxo to +cxo. We 
will replace G{P', P"; r/i, . . . , rjn) by G{P', P"; r/i, . . . , r]*), but maintain the 
dependence on rjn of e*^^", since it belongs to the Fourier representation of 
the delta function, a non-continuous and singular function indeed. We have 
two stationary points ai rjn = and rjn = —2w/{3z) = — (-P' + P")/2 These 
two values are indicated as ai and ^2 respectively. The corresponding values 
of the second derivative are h"{ai) = 2w and h"{a2) = —2w. Using the 
procedure illustrated above, we obtains: 

G{P',P";r]i,r]2, . . . , ai)[227r/t/i"(ai)]i/2e**^("i)+^^°i-^^'/(2*'^"("i)) + 
G{P',P";r]i,r]2, . . . , a2)[2ivr/t/i"(a2)]^/^e^*'^(°2)+ifca2-»fcV(2t/^"{a2)) 

Integrating now on rjn-i the two terms of the expression above we obtain the 
sum of four terms, which are again doubled after the next integration and so 
on. Lastly the asymptotic expression of Ip/p// become 
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In the multiple sum, each oiii, i2, . . . ,in assumes the value 1 or 2, correspond- 
ing to the two values of the stationary points for each relative momentum. 
Finally, performing the integration on k, we obtain: 






-^)-^ E'G(^'^ ^";«n, «.„..., aJ-S_====n 
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^iTJ^im^) i = 1 ^-M'{a, 



where Yl! means that the terms for which I]?=i l/^"('^jj) = ^-^e excluded. 
Such terms indeed give zero after the integration on k. After substitution of 



'iY7j=i l/^"('^ij)| with the lower possible value a/|1/(2w)|, we have asymp- 
totically 

ili2---in 

For the primed sum we may also write: 

T.U... IGIP', P"; «,„...,«,„) I < 
2"[Em,....„ |$(P',«n, • • . ,«^J||$*(P",«n, • • • ,«^JI]/P 

where p is the number of terms which contribute with non zero value to the 
sum. We used the definition of G{P', P"; rji, 772, • • • , Vn) already given and the 
fact that p is equal or less than the 2" possible combinations of the indexes 
ij. 

At this point we will compare the off-diagonal matrix elements with the 
diagonal ones. Since the last are constant of motion, we do not expect it 
would change the result, but to obtain a dimensionless expression. To this 
end we need to establish a further bound. Looking at the definition of the 
reduced matrix elements we may state that 



\Pf 



Ipp = J $(P, r^i, . . . , r/„)$*(P, rii,..., 7]^) dS 



is evidently greater or equal to the result of the integration over a finite 
volume as the integrand is a positive quantity. Let —Mc < i]i < Mc be a 
finite integration interval for each variable. Recalling that only n — 1 of them 
are independent, as ?7„ = — Z]j=i Vj'- 



where ui,U2--- denotes an internal point of the hypercube chosen as in- 
tegration domain. We recall that the function \^{P,ui,U2 ■ ■ ■)\ is assumed 
as bounded everywhere. Comparing now the off-diagonal elements to the 
diagonal ones we will consider, to be democratic, the ratio of |Pp,p„| to 
i\Pp/pMPpnp„\Y^'^- Using the inequalities already found: 

Kp,M\{Pp.p,\\Pp.p.\y^'<2A\t/rr'^-~'y' . 

The constant A denotes the quantity: 

[EU....„ MP', an, . . .,a^„)\mP",a^,,. . ■,a^„)\]/p 

mp>,u[,u',...)\mp",uiu'i..)\ 

and r = ^,pi^^^„2< is a characteristic time scale of the process such that when 
t ^ T the absolute values of the off diagonal elements are vanishing. It 
would be pointed out that actually we have decoherence in the basis of (non- 
relativistic) kinetic energy E = -^ of the cm.. Introducing E as a, more 
appropriate variable the decoherence time is better written as r = 3^^^^^™^//! - 
It is inversely proportional to the difference \E' — E"\ and also decreases with 
increasing total mass M as one may expect intuitively. With \E' — E"\ of 
the order of a few eV and with, say, ten particles, the characteristic time is 
very small, about 10^^^ sec. 

As we anticipated and contrary to our initial prediction, the result we 
have just found shows that the reduced density matrix diagonalizes on the 
basis of P^ and not on that of total momentum P. That is, they are the 
matrix elements for which \P'\ 7^ \P"\, which we found go to zero as time 
increases. But our expression for the asymptotic behaviour is not valid if 
P' = —P". In such a case, in fact, the function h{r]) depends only on 77^, 
we have only a stationary point of the second order and the asymptotic 
evaluation of our integral has to be worked out separately. To this end, 
the integration interval in r]i is split into (0, 00) and (—00, 0) and each of 
the two integrals is integrated by parts a number of times, differentiating 
G{P', P"; r^i, ?72, . . . , rjn) and integrating the remaining factor of the integrand. 
Taking only the dominant term into account, asymptotically we have: 

p,pMPp.p„\y^'<B\t/Tr'^-''y'\ / dpiAm-i . 



The constant B is: 

|$(P^O,0,0,...,0)||$(P^O,0,0,0,.■.,0)| 

mp',u[,u',..:)\mp",u'i,u'i...)\ 

and r' = 3„2|^/_^p//|c - For values of |P'— P"| of a few eV/c and always assuming 
n ~ 10, still t' ~ lO^^^sec. The function Ai{f3) indicates the Airy function 
- J(^°° cos{Py + y^/3)dy. We assumed time to be positive and P' > P". The 
asymptotic estimate of the other half of the off-diagonal matrix elements for 
which P' < P" is immediately obtained since \Pp„p, \ = \Pp,p„ \ by hermiticity. 
Ai{P) for large and positive values of (3 tends exponentially to zero. It is 
also going to zero for (3 which tends towards — oo, but now the behaviour 
is like |/5|~^'^. However, for a sufficient number of particles (five or more) 
[Ai{P)]"- is integrable. Summing up, \Pp,p„\/{\Pp,p,\\Pp„p„\Y^'^, if the number 
of particle is as indicated, also tends asymptotically to zero when P' = —P" 
and we have a complete diagonalization on P basis. 

These findings are strongly dependent on the assumption that the system 
of particles is isolated. It presumes a not very realistic situation since from 
physical viewpoint it is hard to exclude external infiuences, and even under 
the slight perturbation the time behaviour we obtained for the reduced matrix 
elements would not persist. So, even if we intended to present our findings 
strictly as a mathematical result, a brief discussion about the persistence of 
the present decoherence effect in the case of external perturbations imposes. 
We will discuss this point in comparison with the case of a free non-relativistic 
particle. 

The time dependence in the momentum basis of a non-relativistic free 
particle density matrix elements is simply: 



i){k')i)*{k")e 
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where ip{k) is the initial state. Under a time average the expression above, 
compared to the diagonal elements, which are constants of motion, becomes: 

with, denoting as T^ the kinetic energy, a "decoherence time" r, ^^ = .^ _^^ , . 
Since the integral gives a limited oscillating term the asymptotic behaviour 



for the correlations of two different eigenstates of kinetic energie is as the 
inverse power of time, namely as {t/T,^J~^. Let now consider a constant 
perturbation given as a generic function of the energy ef{T), which commutes 
with the free Hamiltonian, and defines T^" = T^/ — A. 

Taking the time average requires to substitute the former integral with 
the new expression: 



t Jo 



For each value of A for which the argument of the exponential is zero the 
integral gives just one. Then for each k' there are so many values of k" as 
the zeroes of the function f{Tk') — f{Tk' — A) for which the averaged matrix 
elements do not go to zero with increasing time. In the case of decoherence 
we considered the role of a free particle is played by the centre of mass. 
Adding any constant perturbation ef{Tp) to the Hamiltonian of our isolated 
system will not disturb at all the asymptotic behaviour of any reduced matrix 
element. In fact we merely have a change of their phases, whereas their 
absolute values are left unchanged. This is a rather artificial consideration, 
since the sought disturbance is not very realistic. However this example may 
deserve to show that our case of decoherence is persistent under disturbances 
which only touches the matrix elements phases. 

Coming back to the case of the free non-relativistic particle, we now con- 
sider a perturbation due to a constant force, adding to the free Hamiltonian 
the term H' = xF, which may be written as xmc/ri. Choosing the acting 
force as the weight force, the quantity ti is about 3 x lO^sec. 

It is easy to see that the wave function of the free falling particle evolves 
with time in momentum basis as: 

'4'{k + mct/ri) e ™' ^' ^^i 

with the same initial condition as before. The above expression means that 
in addition to the phase variation with time of each plane wave component, 
the wave packet "moves" in momentum space at the rate of itlc/ti per sec. 
This entails that, for a square integrable wave function, ilj{k + mct/Ti) goes 
in the long run to zero for any fixed value of k and hence any density matrix 
element with fixed indexes do the same. For a wave packet with initial spread 
Ak equal to some fractions of mc it will happen after a time of the order of 
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year. So we may disregard this kind of time dependence for the first minutes. 
Within this approximation the time average of the normahzed density matrix 
elements gives: 

|V^(A;0||^(A;")| t Jo 
We rewrite the time factor as 



,2 



— / e ■''■'='= "1 at 
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where r, = .„ ^^ , was already encountered and r^ = —^ — ^(r, nY^'^. 

free ^•y—'J'y,\ -J '"■ mc ^ tree <- 1 

We point out that t^ is generally greater than r^^^^. The asymptotic expan- 
sion gives: 



*2 , —itlr, „i ,2 
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The first dominating term is easy written with a couple of Fresnel integrals 



C{z) and S{z) where z = J2/Trt/Tm- 

/ e ^rft=^--— [C(v/2At/rJ-^5(v/2At/rJ] 



t JO t/r„ 

For z < 1 the Fresnel integrals are increasing functions of the argument, and 
both of them goes to the constant value 1/2 in the limit of large z. In conclu- 
sion in the case of non-relativistic particle subjected to weight force we have 
still decoherence with a 1/t asymptotic behaviour. The decoherence time 
however is now the greater quantity r^ instead of t^.^^^. We also examined 
our system of free quasi- relativistic particles, in the case \P'\ ^ \P"\, under 
the perturbation of the weight force applied to the cm.. We encountered 
again three characteristic time scales, one of which is the already defined 
r. At time t ^ t but shorter than the analogue of r^ and in the worse 
situation the off-diagonal density matrix elements differ from zero. We have 
still decoherence, but with a greater decoherence time, the analogous of r^. 
However the time behaviour for the decays of correlations is now as 1/t"^^. 
We may conclude that the decoherence mechanism we presented here is 
as robust or a little more robust to perturbations than the non-relativistic 
free particle case, i. e. not very robust. 

11 



For the sake of completeness, we will also consider the case of relativis- 
tic particles, using the so-called "relativistic Schroedinger equation". This 
equation has been used in various situations, for example to study the sta- 
bility of matter constituted of fixed nuclei and electrons and the collapse of 



stars of self-gravitating bosons or fermions||20||. It is obtained replacing in 
the usual Schroedinger equation the kinetic energy of a particle with its full 
relativistic expression. Let us suppose that there is a Lorentz frame such that 



{"Hi) p ^ l-P/''^! fo^ ^ set R of values of P, the average value in question being 
calculated by means of the above mentioned function <I'(P, r/i, 772, . . . , rjn)- In 
these conditions we can use a procedure similar to the previous one to calcu- 
late the time behaviour of the off-diagonal reduced density matrix elements 
whose indexes P', P" refer to eigenvalues of the total impulse belonging to 



the set R. In fact, the kinetic energy Jrn?c'^ + {P/n + r/j)^ of each particle i 
can be developed with good approximation in power of r/j retaining the terms 
up to rff. Repeating, mutatis mutanda, the previous steps we once again find 
that the off-diagonal elements of the reduced density matrix compared to 
diagonal ones decrease asymptotically with time, as follows: 

The constant C is: 

(|$(P^^;i,^;^..)||$(P",<,^;^..)|)V2 

denoting as 71,72 the two new stationary points. The decoherence time is 
now 



nM^c^\E" - E"\{E"'^ + E'E" + E"'^) 

where E'^ = M'^c^ + P^d^. For E' ~ E" ~ lOMc^, \E' - E"\ ~ leV and 
considering about ten particles, we obtain r^ ~ 10"^'^ sec. All of this applies, 
we remind, if both \P'\ and \P"\ belong to R. As before, the result is that the 
reduced matrix diagonalizes on the basis of P^. But here too our expression 
for the asymptotic behaviour is not valid if |P'| = \P"\. In that case we have 
a double stationary point, the evaluation of the off diagonal reduced matrix 
elements has to be worked out separately with the result that asymptotically 
\Ppip„\ depends on time as \t\^^"^^^^^. 
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4. CONCLUDING REMARKS 

Using non-relativistic quantum mechanics, and within a reference frame 
in which the system is on average at rest (i.e. its average total momentum is 
zero), we have shown that a one- dimensional system of a sufficient number 
of free quasi-relativistic particles, even if isolated, may decohere on the basis 
of the total momentum if we considers correction to kinetic energy to terms 
of second order in 1/c. In other words, we have established, under the con- 
ditions stated above, a superselection rule for the total momentum. For the 
sake of completeness, we obtained similar results using the full relativistic 
expression for the kinetic operator for those Lorentz frames in which the rel- 
ative momenta are much smaller than P/n, where P is the total momentum. 
The decoherence we are talking about was obtained tracing upon the rela- 
tive momenta. Our choice of that coarse-graining seems natural if we want 
to look at the system as a whole. The higher the number of particles, the 
faster the decoherence, meaning that it becomes more effective in the macro- 
scopic limit. What we found in the case of relativistic corrections within a 
non-relativistic quantum mechanics framework may reasonably suggest that 
isolated systems of charged or gravitational interacting particles, for which 
the approximate Hamiltonian we mentioned is adopted, present decoherence 
on the basis of total momentum if we look at the system as a whole and 
ignore the internal degrees of freedom. Systems of this sort might include 
atoms with many electrons, for instance. True, the Hamiltonian in ques- 
tion is not only approximated but also ignores spin, so it may be considered 
"semi-classic" . The presumed effect is interesting anyway. If the number of 
particles is very big, such that the system may be considered macroscopic, 
it cannot be supposed to be isolated, due to the extraordinary density of 
its energy levels and to the action of surrounding objects, even if this tends 
towards zero. Also we learned that our decoherence mechanism is not bery 
robust to perturbations and so other non relativistic interaction terms will 
cause a more effective and different phenomenologically decoherence of the 
cm. variables. The action of the external environment will certainly prevail, 
so that the reduced density matrix for the collective variable of our system 
will diagonalize on the basis of the position, as shown in many examples [11 



But a mesoscopic object, with a very small number of components, may 
probably be considered isolated, if only for a short period of time. In which 
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case we get the following situation: macroscopic bodies appear in a localized 
state, while particles made up of a small number of components would ap- 
pear in well defined states of the total momentum. Further, while not being 
realistic, the model we have discussed for fully relativistic particles can be 
applied to the case of multiple production in the two bodies collision at high 
energy. As it is well known, the secondary particles produced in these exper- 
iments have a transverse momentum pt, which is negligible compared with 
the momentum parallel to the collision axis (pi), and in fact, one-dimensional 
variables have been introduced already to describe the processl^TJ. A very 



popular phenomenological model considers the secondaries in the final state 
as non-interacting particles grouped in clusters having a classical statistical 
distribution in the longitudinal momentum or rapidity [^. Let us think of a 
cluster in terms of our model (adapted for the purpose) and consider the rel- 
ative momenta of the secondaries emitted from it to have values ~ pt, which 
is much less (outside the central zone) than the total momentum of the sys- 
tem. Then we must conclude that the reduced density matrix of clusters 
rapidly diagonalizes on the basis of their total momentum. Which justifies 
their description in terms of classical statistical distribution. 
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